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ABSTRACT 

We compute spectral energy distributions (SEDs) for passive T Tauri disks 
viewed at arbitrary inclinations. Semi-analytic models of disks in radiative and 
hydrostatic equilibrium are employed. Over viewing angles for which the flared 
disk does not occult the central star, the SED varies negligibly with inclination. 
For such aspects, the SED shortward of ~80 fim is particularly insensitive to 
orientation, since short wavelength disk emission is dominated by superheated 
surface layers which are optically thin. The SED of a nearly edge-on disk is 
that of a class I source. The outer disk occults inner disk regions, and emission 
shortward of ~30 fim is dramatically extinguished. Spectral features from 
dust grains may appear in absorption. However, millimeter wavelength fluxes 
decrease by at most a factor of 2 from face-on to edge-on orientations. 

We present illustrative applications of our SED models. The class I source 
04108+2803B is considered a T Tauri star hidden from view by an inclined 
circumstellar disk. Fits to its observed SED yield model-dependent values 
for the disk mass of ~0.015 M and a disk inclination of ~65° relative to 
face-on. The class II source GM Aur represents a T Tauri star unobscured by 
its circumstellar disk. Fitted parameters include a disk mass of ~0.050 M Q and 
an inclination of ~60°. 



Subject headings: circumstellar matter — radiative transfer — stars: pre- 
main sequence — accretion, accretion disks — stars: individual (GM Aur, 
04108+2803B) 



1. INTRODUCTION 



Excess infrared (IR) emission from T Tauri stars is thought to originate from 
circumstellar disks (Mendoza 1968; Shu, Adams, & Lizano 1987; and references therein). 
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Passive disks are the simplest to consider. By definition, they lack intrinsic luminosity and 
reradiate the energy they absorb from the central star. 

Hydrostatic, radiative equilibrium models for passive T Tauri disks are derived by 
Chiang & Goldreich (1997, hereafter CG). The disk surface flares outward with increasing 
radius and intercepts more stellar radiation, especially at large distances from the central 
star, than a flat disk would (Kenyon & Hartmann 1987). An optically thin layer of 
superheated dust grains ensheaths the entire disk. Dust grains in this surface layer are 
directly exposed to central starlight and reradiate to space about half the stellar energy 
they absorb. The other half is emitted towards the midplane and regulates the temperature 
of the cooler disk interior. Vertical temperature gradients in externally illuminated disk 
atmospheres are calculated in detail by Calvet et al. (1991), Malbet & Bertout (1991), and 
D'Alessio et al. (1998). 

CG compute spectral energy distributions (SEDs) of passive disks viewed face-on. 
The calculated SED is fairly constant over the thermal IR, in accord with the observed 
flatfish excesses of T Tauri stars such as GM Aur. Spectral features from dust grains in the 
superheated layer appear in emission when the disk is viewed face-on, similar to solid-state 
emission lines evinced by T Tauri stars and Herbig Ae/Be stars (Cohen & Witteborn 1985; 
Waelkens et al. 1996; Waters & Waelkens 1998). 

What is the spectrum of a passive disk whose midplane is inclined at an arbitrary 
angle, 9, to the plane of the sky? In §|2], model assumptions and results pertaining to the 
continuum SED are set forth. In §|3|, we study how spectral signatures of some dust grain 
resonances change from emission to absorption as the disk is viewed increasingly edge-on. 
Applications to observations are contained in §[|. There, we assess the possibility that the 
relatively low 12 fim flux of GM Aur might be caused by a non-zero disk inclination, as 
opposed to a central AU-sized gap. We close by considering whether differences between 
class I and class II SEDs might reflect differences in viewing angle rather than evolutionary 
status.0 



2. INCLINATION DEPENDENCE OF CONTINUUM SED 



1 In the classification scheme of Lada & Wilking (1984; see also Lada 1987), class I sources exhibit SEDs 
which rise from 2 to 10 /xm, i.e., —3 <J n <J where vF v oc v n . For class II sources, <J n <; 2. 
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2.1. Model Assumptions 

We consider a passive disk in radiative and hydrostatic balance around a T Tauri star. 
The model employed is identical to that derived in §2.1 and §2.3 of CG. Symbols and values 
of free parameters are listed here in Table 1, and henceforth are used without explanation. 
We measure position by cylindrical radius, a, and vertical distance above the disk midplane, 
J 

Derived disk properties relevant to our present investigation are summarized as follows. 
At each radius we distinguish two regions: the superthermal surface layer which is directly 
exposed to light from the central star, and the cooler, diffusively heated interior which the 
surface encases. Variables evaluated in the former region are denoted by a subscript s, 
while those in the latter region carry a subscript i. Pertinent results from CG include the 
dust temperature at the surface 

r**^K, (i) 



2/5 
a AU 



and the gas density in the interior 



i £ ( z 2 \ 
Pai = ~~fKz. T ex P I — oli • ( 2 ) 



Both the gas scale height, h, and interior temperature, Tj, take on different forms depending 
on the vertical optical depth of the interior. Inside 84 AU, the interior is opaque to its own 
reprocessed radiation, and 

Ti ~ ~577 K , (3a) 



3/7 
a AU 

h/a « 0.04 Oau- (3b) 

Between 84 and 209 AU, the interior is optically thin to its own radiation but still thick to 
radiation from the surface; here 

Ti w 21 K , (3c) 



2 The disk radius in astronomical units is denoted by q,au- 
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Table 1. Free Parameters 


Symbol 


Meaning 


Value 


M* 


Stellar Mass 


0.5M Q 


R* 


Stellar Radius 


2.5R Q 




Stellar Effective Temperature 


4000 K 


£ 


Disk Surface Density 


10 3 aXu /2 gem" 2 


a. 


Inner Disk Radius 


6R* = 0.07 AU 




Outer Disk Radius 


2300072* = 270 AU 




Dust Opacity at Visual Wavelengths a 


400 cm 2 g _1 




Grain Emissivity 


f 1 if A < 2nr 
\ (Iff otherwise 


r 


Grain Radius 


0.1 /xm 


P 


Grain Emissivity Index 


1 



a Absorption by dust grains is assumed to dominate the continuum opacity 
from visible through millimeter wavelengths. 
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h/a ~ 0.15 V2 . (3d) 

Finally, in the outermost regions of the disk, the encased material is transparent to radiation 
from the surface, and 



/ 209 \ 19/45 

7^21 K, (3e) 

Va A u/ 



OAD \ 1V45 



V-0.23 1^) . (3f) 



Boundary conditions are as follows. The superheated layer is located at \z\ = H « Ah. 
Its visual optical depth normal to its flared surface is equal to the angle a at which stellar 
rays penetrate the disk (refer to Figure 3 of CG): 



OAR, d (H\ 

+ - ■ ( 4 ) 



da V a / 

Radially, the disk extends from the silicate condensation boundary at a i: to an outer radius, 
a a , at which H w a. 

The occulting angle, 



= arctan — - = 45° , (5) 

is a natural angle with which to compare the viewing angle 9. For 9 <^ 9q, both the central 
star and most of the disk surface are in direct view. For 9 ^ 6 , star and inner disk are 
blocked from view by the flared outer "wall" . 



2.2. Interior SED vs. 

For clarity, we first consider only the disk interior. Figure [l] depicts the viewing 
geometry and coordinate system. Define x and y to be orthogonal spatial axes in the plane 
of the sky, centered on the star. The SED is the integral of specific intensity over projected 
disk area, viz., 

ra ry{x) 

L v = 4ird vF v = 8iru / dx I dy I v , (6) 

Jo J~y{x) 




Fig. 1. — Schematic of viewing geometry and coordinate system. The disk interior is capped 
above and below by the superheated surface. To calculate the SED, specific intensities along 
lines-of-sight (dotted lines parallel to I) are summed over the projected disk area in x and y. 
Here 6 is slightly less than the occulting angle 9 . 
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where 



2/0) = \f> 



a 2 — x 2 cos 9 + H(a ) sin # 



(7) 



traces the disk boundary on the sky, d measures distance to the source from Earth, and 
I u (x,y) is the specific intensity; 



Here B v is the Planck function, and /, / both measure line-of-sight distance from the 
observer. Since temperature, density, and opacity are specified by our model, calculating 
the SED at arbitrary inclination is primarily an exercise in multidimensional integration. 
We perform these integrals numerically using standard techniques (Press et al. 1992). 
Optical depths are evaluated using a Romberg integrator with a mandated fractional 
accuracy of 10 -3 . Specific intensities are computed on an adaptive stepsize grid, where the 
source function (B u ) changes by no more than 20% between adjacent grid points. Integrals 
over y are performed using either a Romberg integrator or, in cases where the integrand is 
very sharply peaked and the Romberg integrator takes prohibitively many steps, a 60-point 
Gaussian quadrature routine. The final integration over x employs the trapezoidal rule on 
a 60-point logarithmic grid. Answers are routinely checked for robustness against changes 
in stepsize. We verify that the SED for 9 = 0° computed using the multi-dimensional 
integration code matches the SED for face-on disks computed using the 1-dimensional 
integrator of CG. 

Figure ^| displays a family of interior SEDs for different 9. First consider SEDs for 
angles 9 < 9 (9 = 0°,30°,45°). Radiation at wavelengths shorter than ~100 /mi comes 
mainly from optically thick regions; consequently, L v is approximately proportional to the 
areal projection factor cos 9. At longer wavelengths, the radiation is emitted by increasingly 
transparent material, so that L v tends to be independent of orientation. 

For angles 9 > 9q{9 = 65°, 85°), short wavelength emission from the inner disk is 
strongly extincted by the intervening outer disk. However, millimeter (mm) wavelength 
fluxes only drop by a factor of 1.8 as the disk is viewed increasingly edge-on, since about 
half of the outer disk remains visible in that limit. 





with 




(9) 
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Fig. 2. — SED of the disk interior only, as a function of viewing angle. Short wavelength 
emission from the inner disk is extinguished by the intervening outer disk for 9 > 9 . 
denotes the luminosity of the central star. 
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2.3. Superthermal + Interior SED vs. 

Next we compute the total SED. The superheated layer is treated as a plane-parallel 
atmosphere having visual optical depth a normal to the surface. Whenever a line-of-sight 
intersects the surface, we increment the specific intensity from the interior byi 



AI U = B U {T C 



ds) 



as. 



exp 



\n ■ l\ 



exp(-r^) . (10) 



Here, n and I are unit vectors normal to the surface and parallel to the line-of-sight, 
respectively, e s is the Planck-averaged dust emissivity at the surface, and t Vj i is the 
intervening optical depth between disk surface and observer. 

Figure displays a family of total disk SEDs labeled by 9. The behavior of the SED 



for A ^ 80 /im is similar to that described in § |2.2| , since the interior dominates emission at 
those wavelengths. 

Shortward of ~80 fim, emission from the superheated surface makes a qualitative 
difference to the appearance of the SEDs. One difference is that for angles 9 < 8 , the 
flux is remarkably insensitive to viewing geometryE This is because shorter wavelength 
radiation emerges primarily from the superheated surface layers, which are optically thin 
along most lines of sight. 

For angles 6 > 6q, the intervening outer disk dramatically extinguishes surface radiation 
at the shortest wavelengths. However, the total disk SED between 30 and 80 /im is relatively 
robust to changes in inclination; it decreases by less than a factor of 3 between face-on and 
edge-on orientations. Radiation at these wavelengths originates from surface layers at large 
radii and is only slightly obscured by the tenuous interior at still larger radii. 

Finally, in Figure |], we add the contribution from the central star to the SED. The 
stellar flux fades rapidly with increasing 9 beyond ~55°. 



3 This prescription fails for lines-of-sight which graze the surface tangentially. However, their fractional 
contribution to the integrated flux is of order y^h s /R, where h s « (<91n p g i/dz)~ 1 \ z =H is the geometrical 
thickness of the surface, and R is the radius of curvature of H . Since this fraction is less than 10%, for 
computational simplicity we let the intersected optical depth of the surface be ae s x min(|n • 20). 

4 Slight reductions in mid-infrared flux still occur for 9 < 6$, both because surface layers intermediate in 
radius are seen through the outer disk's atmosphere and because emission from the optically thick interior 
is proportional to the projected area. 
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Fig. 3. — SED of the entire disk (interior + surface), as a function of viewing angle. For 
9 < 9 , the SED is insensitive to viewing aspect since short wavelength radiation originates 
from the optically thin surface, while long wavelength radiation is emitted by the optically 
thin interior. For 9 > 9 , the intervening outer disk dramatically extinguishes surface 
radiation at the shortest wavelengths. 
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Fig. 4. — SED of the entire disk + central star, as a function of viewing angle. At 9 = 60°, 
the visual optical depth to the star is r^y = 4.5. At 9 = 65°, r*y — 22. At 9 = 85°, 
T*,y = 7.6 x 10 5 . 
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3. Spectral Dust Features 

When viewed nearly face-on, passive disks exhibit emission lines associated with dust 
grain resonances. These features arise from material along the line-of-sight which is optically 
thin: line-to-continuum flux ratios are of the same order as the percentage enhancement in 
line opacity. Emission lines may originate from both the optically thin superheated surface 
(as discussed in CG) and the optically thin interior. 

How do such lines vary as the disk is viewed increasingly edge-on? Might they appear 
in absorption instead, as Cohen & Witteborn (1985) find for 7 of the 32 T Tauri stars they 
survey spectrophotometrically? To address these questions, we add 6 Gaussian-shaped 
spectral resonances to our grain emissivity law. In choosing resonant wavelengths and 
strengths, we aim to illustrate the range of effects which occur with variable viewing angle. 
Roughly, our 6 lines can be associated with crystalline H2O ice at 60 and 100 /iin (Waters & 
Waelkens 1998), amorphous silicates at 10 and 20 /im (Mathis 1990), and small carbon-rich 
grains or PAHs (polycyclic aromatic hydrocarbons) at 3.3 and 6.2 /zm (Waelkens et al. 
1996). 

Solid-state line spectra at various disk inclinations are displayed in Figure |5|. Emission 
at all wavelengths varies little with viewing geometry for angles 9 <; 9 . This is true 
particularly at the shortest wavelengths where lines arise largely from the optically thin 
surface layers. At wavelengths progressively longer than ~30 /im, the optically thin interior 
at large radius contributes increasingly to line emission. The drop in the long wavelength 
continuum with increasing inclination reflects the cos 9 dependence of radiation emitted by 
the optically thick interior. 

Emission lines at ~60 and 100 /im persist at angles 9 ^ 9q, although their fluxes are 
lower by a factor ~2 for edge-on as compared to face-on disks. For edge-on orientations, 
only about half of the optically thin interior remains visible. Emission around 20 fim is 
more robust to changes in inclination; most of this radiation emerges from the superheated 
surface at large radii and passes through the rarefied disk interior at yet larger radii with 
little attenuation. 

At the shortest wavelengths, where occultation of the inner disk is significant, lines 
indeed appear in absorption, though only over a limited range of inclinations between 9q 
and 90°. For a dust grain resonance to appear in absorption, disk inclinations must be 
sufficiently high that line emission from the superheated surface is extinguished. However, 
9 cannot be so high that line and continuum optical depths both exceed unity at large 
radius where the disk interior is nearly isothermal. Only over an intermediate range of 
inclinations does radiation on and off the resonant wavelength probe a variety of interior 
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Fig. 5. — Solid-state line spectra at various disk inclinations. Resonances at the shortest 
wavelengths where opacities are highest appear in absorption over a limited range of viewing 
angles between 9 and 90°. Longer wavelength features persist in emission. 
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disk temperatures. For the 3.3 /im resonance to appear in absorption, 60° <^ i <^ 65°. At 
6.2 fim, opacities are lower, and inclinations needed for absorption higher: 65° <^ i <^ 70°. In 
principle, near-infrared line spectra may provide a sensitive diagnostic of disk inclination. 

While the trends outlined above apply qualitatively to all passive disks, quantitative 
conclusions are model-dependent. For example, in our standard model, no absorption line 
at 10 fim appears because line-of-sight optical depths are too small to extinguish emission 
from the superthermal surface. However, such an absorption line would appear in the SED 
of a sufficiently massive, inclined disk. 

4. Applications to Observations 

4.1. GM Aurigae and Central Holes 

GM Aur is a particularly clean, well-studied T Tauri system to which we can apply 
our SED models. The single central star has mass O.72M , luminosity O.7L , and age 
2 x 10 6 yr (Beckwith et al. 1990). Aperture synthesis maps in 13 CO(2 — > 1) evince a 
rotating circumstellar gas disk having approximate projected dimensions 950 x 530 AU, 
inclined at 9 ~ 30° (Koerner, Sargent, & Beckwith 1993). Hubble Space Telescope (HST) 
images in scattered visible light reveal the disk surface to be flared, and suggest % is closer 
to ~60° (Koerner 1997; Stapelfeldt et al. 1995). These same observations detect no outflow 
or remnant envelope surrounding this relatively evolved system. 

CG model the SED of GM Aur, but with a face-on disk. They suggest that a non-zero 
disk inclination might cause line-of-sight radiation from the inner disk to be absorbed by 
the outer disk, thereby explaining the relatively low 12 /xm flux without invoking a central 
hole. Here we investigate this possibility by refitting the SED with an inclined disk. 

Figure ^| displays our new fit with a disk inclined at 9 = 59°. Despite the non-zero 
inclination, none of our fit parameters changes significantly from those given by CG for the 
face-on case. Fluxes at IRAS (Infrared Astronomical Satellite) wavelengths emerge mainly 
from the optically thin surface layers, and those at mm wavelengths emerge mostly from 
the optically thin interior; neither is sensitive to viewing angles 9 near or less than 9q. ^ 

Our model reproduces the observed modest extinction to the central star of Ay = 0.5 
magnitudes, but at the cost of still requiring a large central gap having a radius ~60 times 
greater than the dust sublimation radius: cij = 4.8 ± 2.7AU, where the uncertainty reflects 



5 9o for this disk equals 51°. 
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Fig. 6. — Observed SED (filled circles) of GM Aur and accompanying inclined passive disk 
model. Fit parameters are as follows: 9 = 59°, £ = 2x 10 3 a~^j 2 g cm -2 , (3 — 1.4, r = 0.3 /im, 
a Q = 311 AU (which implies 9q = 51°), and aj = 4.8 AU. The formal requirement of a large 
AU-sized central gap persists despite accounting for a non-zero viewing angle. 
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that of the 12 /im IRAS point. There is no inclination for which the outer edge of our 
model disk attenuates the 12 /im flux without also extinguishing central starlight by several 
magnitudes. Line-of-sight column densities to surface regions inside ~5 AU (where most 
of the 12 fim emission originates) are only ~20% greater than the column density to the 
star. This small difference in the amount of obscuring material is completely insignificant; 
the opacity at 12 /im is 13 times lower than at visible wavelengths. Thus, the low visual 
extinction to the central star implies that we should have a clear view of disk regions inside 
a few AU at near-infrared wavelengths. 

In spite of these considerations, we remain skeptical of the existence of such a large 
central gap devoid of dust. A more palatable alternative might be that the inner disk's 
aspect ratio does not increase monotonically with radius. Such undulations in the height of 
the surface would shadow annular regions from the central star and lower their temperatures. 
In addition, those ripples nearest the observer would hide their warmer, starlit sides when 
viewed at non-zero inclination. Accretional heating may be responsible for such changes in 
surface geometry, as mentioned in CG (see their §3.5); as estimated there, enhancements 
in disk thickness due to heating of the midplane may become significant inside a few AU. 
For an accretion disk which derives its luminosity solely from local viscous dissipation, 
the aspect ratio indeed decreases with radius whenever opacities increase steeply with 
temperature, as demonstrated by Bell et al. (1997). 



4.2. Class I Sources As Inclined Class II Sources 

Might some class I spectra represent T Tauri stars obscured by inclined disks? The 
possibility should be entertained for sources such as 04108+2803B; Figure demonstrates 
that a passive disk inclined at 65° provides a reasonable fit (to within factors of 2) to this 
class I SED. Fit parameters are similar to those of our standard model, and are listed in the 
figure caption. For this particular model, near-IR emission at J, H, and K is interpreted 
as central starlight extincted by the disk's outer edge; inclusion of starlight scattered off 
the disk surface would imply a larger inclination. To highlight the contribution from the 
superheated surface, we also plot the SED with the surface emission removed. Fluxes 
between 3 and 60 /im arise primarily from the superthermal surface, seen both through and 
over the lip of the outer disk atmosphere (refer to Figure [I]). 

Clearly, a class I SED does not imply a unique distribution of circumstellar material. 
High angular resolution images provide additional clues. Many class I sources exhibit 
near-IR nebulosity on scales ranging from 1500 to 3000 AU (Kenyon et al. 1993; and 
references therein); for these sources, an additional, non-disk component of dust is needed. 
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Fig. 7.— Observed SED (filled circles) of class I source 04108+2803B (Kenyon, Calvet, & 
Hartmann 1993) and accompanying inclined passive disk model. Fit parameters for the disk 
are as follows: 6 = 65°, S = 0.3 x 10 3 a A u /2 g cm" 2 , f3 = 1.2, r = 0.1 ^m, a Q = 270 AU 
(which implies 9 = 46°), and a« = 0.07 AU. Stellar parameters are identical to those of our 
standard model T Tauri star. 
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Kenyon, Calvet, & Hartmann (1993) model class I SEDs using rotationally flattened, 
infalling envelopes that are passively heated by central stars. These envelopes extend out 
to larger radii than do our disks, typically 3000 AU, but contain roughly the same amount 
of mass, about O.1M . Bipolar holes empty of material, presumably evacuated by outflows, 
are invoked so that central starlight may scatter off cavity walls towards Earth, thereby 
explaining the observed near-IR fluxes. In imaging observations of three class I sources 
using HST/NICMOS (Near Infrared Camera and Multi-Object Spectrometer), Padgett et 
al. (1998) discover nearly edge-on, flared circumstellar disks having diameters 300-750 AU. 
Near-IR emission is observed to be scattered not only off upper and lower disk surfaces, but 
also off "dusty material within or on the walls of outflow cavities." 

In general, a combination of an inclined, passively heated disk, and a dusty 
bipolar outflow or partially evacuated envelope may best describe class I sources. Since 
04108+2803B reveals no optical or near-IR emission beyond ~140 AU from its central star 
(Kenyon, Calvet, & Hartmann 1993), it is consistent with being a limiting example of a 
simple inclined disk. Two other examples for which this limiting case scenario may also 
apply include the embedded sources 04295+2251 and 04489+3042. Both sources exhibit 
nearly flat excesses between ~3 and 100 /iin which may arise from flared circumstellar disks. 
Neither source betrays extended near-IR emission on scales greater than a few hundred AU, 
or high- velocity molecular gas from an outflow (Kenyon et al. 1993). 

Near-IR polarimetry also addresses the possible presence of envelopes. Whitney, 
Kenyon, & Gomez (1997) determine that linear polarizations greater than ~20% at J, H, 
and K characterize the largest reflection nebulae associated with class I sources. Such large 
polarizations are interpreted as arising from the scattering of central starlight off cavity 
walls in envelopes. In contrast, 04108+2803B evinces fractional polarizations of 5.1% and 
1.6% at H and K, respectively (Whitney, Kenyon, & Gomez 1997) — levels more comparable 
to those in class II sources than in truly younger, more embedded protostars, and which 
may arise from starlight scattered off the disk surface. Resolved polarimetry maps of 
04108+2803B can test our hypothesis. 

Degeneracies inherent in models deduced from a SED may also be broken by kinematic 
studies. Interferometric radial velocity maps in molecular lines may distinguish between 
an infalling envelope (e.g., Hayashi, Ohashi, & Miyama 1993; Ohashi et al. 1997), and a 
rotating disk (e.g., Koerner & Sargent 1995). 

Finally, we note that the inverse problem to 04108+2803B is presented by the star 
HK Tau B: a system observed to possess an edge-on (i « 85°) circumstellar disk and no 
observable envelope, but whose infrared SED is not well measured because of confusion 
from an infrared-bright companion (Stapelfeldt et al. 1998; Koresko 1998). We await 



-19- 



SIRTF (Space Infrared Telescope Facility) which can provide both the SED between 3 and 
180 microns and also images of the superheated disk surface in thermal emission. 

Financial support for this research was provided by NSF grant 94-14232, NASA grant 
NAG5-7008, and an NSF Graduate Fellowship. We thank an anonymous referee for a 
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